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The application of artificial earth satellites in geodesy can be 
either dynamic or geometric. 

In the dynamic approach, an evaluation of the perturbations in 
the orbital motion of a satellite ean furnish data from which the coef 
ficients of the spherical harmonic expansion of the potential function 
can be determined. Using this approach, valuable information on the 
gravitational field of the earth can be obtained. 

In the geometric treatment, the artificial earth satellite is 
used as a reference point or a triangulation station in the sky. In 
this way the positions of observation stations on the surface of the 
earth can be determined. Some aspects of the geometric approach will 
be discussed in this paper. 

The purpose of this thesis is to serve as a demonstration of the 
material studied during its preparation. I would like to emphasize that 
the reader should not construe the methods and procedures used herein 
as being of a standard type because of the following reasons: 

lo The computed orbit from three observations at a 

given station at short intervals is assumed to be correct 
and accurate3 

20 It is assumed that the observed coordinates of the 

satellite taken from Smithsonian Institution Astrophysical 
Observatory Special Report No. 66 and the computed orbital 
elements are referred to the same mean coordinate system; 

3. The Gaussian Method used to determine the orbit is an 

approximate methods; 

lh. It is assumed that only two approximations are necessary 


t 


in the differential process of computing the geographic latitude 


ii 





and longitude 0 
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1. INTRODUCTION 


Since the early days of history man has shown a keen interest and 
an investigative concern about the earth, its characteristics, and its 
dimensions. Since the third centmry B.C. when Pythagoras stated that 
the earth was of spherical shape, many attempts have been made to 
measure its size and shape and to reduce the physical parameters and 
features on to maps and charts. Of the many existing earth sciences, 
geodesy is the one that primarily contains this field. - 

The inception of geodetic science can be traced back to Eratosthenes 
(circa 276 = 19); B.C.) who Was the first man to measure the size of the 
earth. His method of measuring the length of a meridianal arc by timing 
the travel of a camel caravan between Alexandria and Syene (Aswan) and 
measuring the angle of incidence of the rays of the sun shining into a 
well at these places was crude but effective. The degree of accuracy 
he obtained was remarkable when we consider the tools he had available. 

The science showed little development between this time and the 
seventeenth century. At this time the advent of the telescope, loga= 
rithmic tables and the method of triangulation (developed by the Dutch 
scientist Snellius in 1615) took place and geodesy was reborn. 

The history of geodesy can be divided into periods based on the 
figure of the earth that was used as a reference during the respective 
periods. The studies and experiments of Pythagoras, Aristotle, and 
Eratosthenes brought the "flat", or plane era to an end and introduced 
the spherical era. Newton and Huygens, in the eighteenth century, 
introduced the ellipsoidal period which fin ally phased into the 
geoidal and telluroidal periods of this century. 
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With the launching of the first artificial satellite and the begin= 
ning of the space age in this decade a new era in the fieid of geodesy 
is unfolding. By using the geoid as a reference surface we can determine 
the shape of the earth, but not its size. The orbiting artificial earth 
satellites will not replace the geoid or the telluroid but will provide 
the geodesist with a convenient reference station or “base~line in the 
sky" that can be used to measure both the size and shape of the earth. 

Since the launching of our first artificial satellite (Explorer 
I) on 31 January 1958, the field of space science has expanded tremen= 
dously and with the successful results of the early and present projects 
the limits in this science are beyond imagination. At the present, 
the projects already underway and those planned for the immediate 
future suggest the following three areas of space activity: research, 
application, and exploration. (2; P.5] 

The applications include meteorology, commmications, and geodesyo 
No new principles, yet to be discovered are necessary; _no major tech= 
nological problems stand in our way. 

In the field of meteorology, the satellites can provide, through 
telemetry, information on pressure, temperature, cloud coverage and 
stom patterns and also assist in predicting such mankillers as 
hurricanes, cyclones, etc. The day to day practical use of weather 
satellites is obvious. 

In commmications, satellite contributions are numerous: 
telephony, long-range radio commmications, and international tele=- 
vision. 

This brings us to the third application of space activity: 


geodesy. Geodetic satellites can serve both research and application 
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3 
interests. As to the latter, navigation and surveying come closest to 
the immediate interests of mankind. The geodesist can use the satellite 
in both a dynamic and a geometric approach. Dynamically, the pertur= 
bations in the orbits of the satellites caused by changes in the earthis 
gravitational attraction can be mathematically computed and analyzed 
and thus provide valuable information to assist in determining the 
gravity field of the earth. By observing the changes in the orbits 
of many satellites over a long period of time, information on regional 
gravity anomalies can be obtained. From these data we can compute 
the shape of our equipotential reference surface, the geoid. Further, 
observations of the nodal motion or the rate of change of the plane 
of the orbit of the satellite gives us information on the flattening 
ofthe earth. Therefore, in the dynamic approach we can determine the 
shape of the earth. 

The size of the earth can be realized in the geometric approach. 
In this process the satellite will serve as a reference point ora 
triangulation station in the sky. Thus we can connect triangulation 
nets or stations that are widely separated by oceans and continents 
or are located in inaccessible areas where the conventional methods of 
triangulation would be awkward and cumbersome. If so desired, the 
existing geodetic systems could be reduced to a common world geodetic 
network or datum. 

In this thesis a phase of the geometric use of the artificial 
earth satellite in geodesy will be discussed. In the geometric 
approach the geodesist can use the satellite in two basic ways. One 
method involves the observation of the satellite simultaneously from 


two or more stations. The other method requires that the orbit of the 





bh 
satellite be precisely known. Then knowing the orbital parameters or 
characteristics and predicting their values for a desired time we 
observe the satellite from any number of stations on the earth. 


The discussion of the problem will be treated in the following 


First, the orbit of the satellite Explorer I (1958 Alpha) will be 
computed from the right ascension and declination of the satellite as 
observed at three different times from the tracking station at Olifants= 
fontein, South Africa. The three observations have been made at very 
short intervals during the same revolution of the satellite. The 


M EFERRED 
values of right ascension and declination are -appamant as ones S555, ites 


TO THE EPOCH 1950.0, REDUCED 

from camera films of Baker-Nunn tracking cameras. The geographic 
coordinates, latitude Co) and longitude { ) ) of the tracking station 
are knom and its height (H) above sea level is given. 

Second, the same satellite will be observed again from a station 
on the island of Curacao in the Netherland West Indies in the Caribbean. 
The geographic coordinates of this station will be assumed to be unknown. 
Then using the computed orbital parameters and observing the satellite's 
right ascension and declination from this station, the geographic co= 
ordinates of this unknown station will be found. 

Third, the forward and back azimuths and the geodesic distance 
between the two stations will be computed. 

This connection between the two points is called a tie and in 


this instance, an inter-continental tie. The distance between the 


stations is approximately 6,000 nautical miles. he-—preeedure—diseussed 
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2o DEFINITIONS AND ELEMENTS OF THE ORBIT 


The motion of an artificial satellite in its orbit about the 
earth follows the same laws that apply to the movement of the planets 
in relation to the sun in the solar system. 

Johannes Kepler (1571 = 1630), a German astronomer and mathema= 
tician, derived the laws of planetary motion from the observations of 
Tycho Brahe (1546 = 1601), Danish astronomer, with whom he worked as 
an assistant. 

Kepler's Laws are: 


I The orbit of each planet is an ellipse with 
the sun at one of its foci. 


II The line (radius vector) joining the sun to 
each planet sweeps over equal areas of its 
ellipse in equal times, 

III The squares of the orbital periods of the 
planets are proportional to the cubes of their 
mean distances from the sun (foci). 

The relationship of an artificial satellite with respect to the 
earth is the same as that existing between the plmets and the sun as 
stated in the above lams. | 

In other words, we can replace the term planet with satellite 
and replace the sun with the earth and the lays are still applicable. 
Now, since we know the basic form of the orbit is an ellipse we can 
easily describe the characteristics of the orbit. These orbital 
characteristics are knowm as parameters and can be considered to be 
constmt quantities for a particular satellite. The first two para~ 
meters are properties of the ellipse itself and these are the semi~ 


major axis, a, and the eccentricity, e. These two terms define any 


ellipse. 








Figure 1. 


Referring to Figure 1., the eteentnieity is the ratio existing 
between the length of the semi-major axis and the distance of the 
foci, c, from the center, o, ors : 

e = C/a 
The relationship e@ = (ae ~ b°)/a?, where b is the minor semi-axis, 
also exists. The eccentricity of a circle is zero and that of a 
parabola is equal to one. The eccentricity of an ellipse is less 
than one but greater than zero. The eccentricity of the orbit of 
a satellite is nearly circular, that is, slightly greater than Zero. 

Two other parameters of the orbit are the angle of inclination, 
i, and the right ascension of the ascending node, Q e These two 
values define the plane of the orbit in space. Referring to Figure 2, 


the angle i is measured from the plane of the equator to the plane of 
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8 
the orbit. It can have a value between 0° and 180°. The angle (2 is 
measured eastward along the equatorial plane from the vernal equinox, 
to the ascending node. It can have any value between 0° and 360°. 

The ascending node is that point where the satellite crosses the equa~ 
torial plane from south to north. 

Another parameter of the orbit is the argument of the perigee, 

w, or the angle measured from the equatorial plane at the ascending 
node to the major semi-axis (the extension of which contains the 
perigee). The angle w is measured in the direction of the satellite's 
motion and gives us the orientation of the orbit in its plane. 

The sixth and last parameter is the time, To. Either the time 
of passage of the satellite through the ascending node or the time 
of passage through perigee may be chosen. The latter time is usually 
the one used in most computations. 

In summarizing: the a and e specify the size and shape of the 
orbits; i and (2 specify the orientation of the orbital plane in space; 
w orients the orbit in its plane; T, furnishes us with the value to 
locate the satellite in its orbital plane relative to the perigee. 

To position the satellite at any specific time we mst resort to 
a system of angular values called anomalies, of which there are three: 


mean M, eccentric E, and true Ve 





Apogee 


\ 


Line’ of 
Apsides 





Auxiliary — 
Circle 


Figure 3. 


Refer to Figure 3. Angle QFN is the true anomaly V, which is the 
angular distance of the satellite measured from the line of apsides 
along the plane of the orbit. Since the orbit of the satellite is 

an ellipse, the radius vector FQ, or the distance of the satellite 
from the foci, will not be constant. Therefore, the angular velocity 
is not constant. This mst be so if Kepler's second law is to be 
fulfilled. If the angular velocity of the satellite was constant, 
the radius vector would sweep at an average rate and the eccentric 
anomaly would equal the mean anomaly. The mean anomaly is the theo- 
retical angular velocity of a satellite moving uniformly through 360° 


of arc per orbital period. 
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10 


The symbols and definitions of terms used in this paper are listed 


be Lowe 


a 


u 


XgVoS ° 


XA, YaAsyZa 


X09YO059%0 


Semi-major axis of the orbit of the satellite 
(Chapters h, 5); 

Semi-major axis of the International Ellipsoid (earth) 
(Chapters 1, 2, 3) 

Semi=minor axis of the orbit of the satellite 
(Chapters lh, 5)3 

Semi=minor axis of the ellipsoid (Chapters 1, 2, 3) 

Distance from center of orbit to foci 

Eccentricity of orbital ellipse (Chapter ),5); 

Eccentricity of International Ellipsoid 
(Chapters 1, 2, 3) 

Flattening of the ellipsoid 

Direction cosines (used together) 

Angle of inclination of the orbital plane 

Gaussian gravitational constant 


Mean daily motion of satellite (velocity of the 
mean anomaly) in revolutions. 


Auxiliary parameter of orbital ellipse (replaces a) 

Geocentric distance of satellite at perigee 
(Chapter 5); 

Time interval: (Tl + 1T3)k (Chapter h) 

Distance of satellite from observer 

Geodesic distance 

Time interval: (T3 = Tj )k 

Auxiliary parameter of orbit 


Geocentric space rectangular coordinates 


Approximate space rectangular coordinates of 
the observer 


Space rectangular coordinates of the observer. 





A,B,C,D,E,F,G 


E 


Various terms used in Sodano's Fourth Method 
Eccentric anomaly of satellite 

Newton's Universal gravitational constant 
Elevation above sea level 

Gaussia determinants 


(j = 1,2,3,455,6) Constants for International 
Ellipsoid used in Sodano's Fourth Method. 


Difference in longitude between two points 
on the ellipsoid 


Mean anomaly of satellite 

Radius of curvature in the prime vertical 
Geocentric distance of satellite 

Time of observation (j =1, 2,3) 
Time of Epoch 

True anomaly of satellite 


Space rectangular coordinates of the satellite. 


Flattening of the ellipsoid (Chapter 6) 
Geocentric Right Ascersion ~ the angular distace 
at the center of the earth measured eastward 
along the equator from the apparent vernal 
equinox to the celestial meridian passing 
through the satellite. 


Apparent Right Ascension - measured at the 
position of the observer. 


Apparent right ascension obtained through 
the approximate coordinates of observere 


Forward Azimuth 

Back Azimuth 

Parametric latitude (reduced latitude) 
Maximum latitude of the geodesic. 


Normal gravity. 
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Geocentric declination ~ the angular distance at 
the center of the earth measured northward or 
southward from the equator to the satellite. A 
satellite below the equator would have a south 
declination and the angle would have a minus valve. 


Apparent declination = measured at the observer's 
position. 


Computed apparent declination obtained through 
the approximate coordinates of the observer. 


Spherical distance 

True spherical distance 

Greenwich hour angle of the apparent vernal equinox. 
Local hour angle of the apparent vernal equinox, 


Geographic longitude (positive westward and 
negative eastward) 


Geocentric distance of observer 


The angular motion of the satellite during the 
interval between observations. Usually expressed 
in radians. 


Geographic latitude (north positive and south 
negative ) 


Geocentric latitude 


Argument of the perigee. The angle measured 
eastward in the orbital plane from the line 
of nodes to the line of apsides. 


Gaussian determinant 


Vernal Equinox = the point of intersection of 
the path of the sun with the equator at which 
the sun is passing from south latitude to north 
latitude. 


Angle at satellite's position measured between 
the line to observer's position and the line to 
the earth's center (parallactic angle at middle 
observation). 


Angle at observer's position measured between 
the line to the earth's center and the line to 
position of the satellite at the time of the 
middle observations 





Right ascension of the ascending node. 
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30 DETERMINATION OF THE MASS OF THE EARTH 
AND GAUSS! GRAVITATIONAL CONSTANT 


The existing relationship between the earth and the satellite 
can be treated as a modified two body problem in which we assume that 
the two bodies are spheres and that their mass distribution is 
homogeneous. Under these circumstances, we can say that they will 
follow the laws of celestial mechanics and gravitation. That is, they 
will attract each other with a force that is proportional to the 
product of their masses and which varies inversely as the square of 
the distances between their centers. Since the mass of the satellite 
cante considered to be negligible compared to that of the earth and 
because of its proximity, it canlte stated that the gravitational attrac 
tion of the earth is the only attracting force influencing the orbit. 
The force acting between the two bodies is: 


F = ke Me Mg 


where r is the distance between the two bodies, me is the mass of ths 
earth, msg the mass of the satellite, and k equals the Gaussian gravita~- 
tional constant. The numerical value of the gravitational constant 
depends on the units of mass, time, and distance chosen. 2, P.57| 
Kepler's third laws 

k* (me + mg) = na? 
is a simple relation between mass, time, and distance that may 
conveniently be used in practical application of elliptic motion to 
find any one of these three quantities when the other two are known. 
The units of mass, time, and distance can be arbitrarily chosen and k 
can either be computed or derived from observations. In applying the 
units to determine k for use in observation of planetary orbits the 


following units are used: the sun's mass as the unit of mass, the 


ly 





15 
ephemeris day as the unit of time, and the semi-major axis of the earth's 
orbit (which approximates very closely, one astronomical mit) as the 
unit of distmce. In the case of the artificial earth satellite a 
revision of these units must be made in order to find a realistic value 
of the Gaussian constant that will be applicable. Bear in mind that 
we are theorizing that the earth is the only body influencing the motion 
of the satellite. This is a sound theory but the other assumptions we 
made = that the earth is a sphere and that the mass distribution is 
homogenecous= are not absolutely correct. In fact, if the earth was 
spherical and the mass distribution was constant throughout, then the 
orbit of the satellite would be nearly perfect and very simple to 
predict. The plane of the orbit would be fixed, the parameters of the 
elliptical orbit would remain constant, and the center of the earth 
would act as one of the foci. 

However, the earth is an ellipsoid of revolution with a bulge of 
mass at the equator and a non—homogeneous distribution of mass through= 
out. These characteristics plus the effect of the earth's atmosphere 
and other perturbations have a marked influence on the behavior of the 
satellite in its orbital travels. 

The effect of the attraction of the equatorial mass of the earth 
is to pull the satellite toward the equator, thus decreasing the inclin- 
ation. Behaving like a gyroscope the actual effect is to cause the 
plane of the orbit to precess westward for an eastward launch about the 
focus. This motion causes the perigee to advance to different points 
along the orbit. The perigee moves in the direction of the satellite 
motion when the inclination is less than 63.°5 and moves inthe opposite 


direction when the inclination is greater than 63.°5. Consequently, 





16 
the elliptic arbit rotates about the focus in the orbital plane. Reduc-~ 
tions of the periods of this perigee motion from the observations provide 
data to determine the shape of the earth. [11, Poli 

The non—homogeneity of the earth's mass distribution causes per~ 
turbations in the orbit of the satellite. By observing a number of 
satellites for long durations of time we can compute the deviations 
from the normal gravity and thus obtain values for the regional gravity 
aWbiielies, 

Proceeding withthe computation of the Gaussian constmt, k, the 
following units will be used: the earth's mass as the unit of mass, 
the ephemeris day (86,100 seconds) as the unit of time, and the semi- 
major axis of the earth as the unit of distance. Compare these with 
the units mentioned earlier in this chapter that are used in conjunction 
with planetary orbits. 

Note that all computations are based on the International Ellipsoid 
(Hayford, 1910) which has an equatorial radius of 6,378,388.000 meters 
and a flattening of 1/297.0. Here again rises an arbitrary value. 

There are many reference ellipsoids from which to choose but the varia~ 
tion in size from one to another is small. The International Fllipsoid 
has been adopted by many countries and it is the one that is used by 
Smithsonian Astrophysical Observatory in their computations, so I 
decided to use this as the reference body. 

The international gravity formula adopted by the International 
Union of Geodesy and Geophysics in 1930 is [35 P.7h] 

f* at B sin? + € sin°2D ) 
where Te ~ gravity value at the equator, D = latitude of dation, 


B = gravitational flattening, € # theoretical constant. The 


Le 
coefficients of the formula are: 
7 = 978.0h9(1 + 0.005288 sin? - 00000059 sin? 2P)c m/sec? 
The gravity value at the equator is the resultant of the gravitational 
attraction and the centrifugal force due to the earth's rotation. 
Therefore , 
Te = w°a = 978.0190 cm/sec? 
where Te = attraction of the earth at the equator and wa = centrifugal 
force. 
If wea = oo )@r where wis the angular velocity of the earth, 


T is the sidereal day in seconds,and r the radius of the ellipsoid then: 


ie )*r = 3.380545 cm/sec“ 


Then 
Te = 3.3805 = 978.090 


Te = 978.0490 + 3.43805 = 981.4295 gal 
where a gal (from Galileo) equals a unit of acceleration equal to 
1 cm/see*. 


The mean density of the earth is computed from the formla 5] : 
( 


l+e é / © 
Po = 271 G C.5as (arctan e ~ eal <a) 


where Po = To/a = attraction of earth at the equator divided by the semi- 
major axis, @G@ = Newton's universal constant of gravitation which equals 
8 


6.673 + 0.003 x 10™ ces, e' = the secondary eccentricity, and Oe 3 
the mean density of the earth. To solve the equation, two of the three 
quantities Te, a, or Cn must be known. In our example, the Te has 
been determined and the value of the semi-majar axis of the Interna= 
tional Ellipsoid is used. The computed mean density is (see Table 1): 
O m.® 5251300437 em/am 
The volume of the International Ellipsoid is computed from the 


formulas 
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v=h/3 7 a@b = 1.083320781 x 10°? cw 
where b = 6,356,911.9),62 meters, the semi-minor axis of the ellipsoid. 
Using these values, the computed mass of the earth is: 
M = VO =(1.083320782 x 20°) (5.51300437)= 5297235220 x 10°’ gm. 
To obtain the proper unit for k in terms of mass of the earth for 
unit mass, time in days, and distance in terms of the earth's semi- 


major axis the following conversion is made: 


i nae com> a: 86,100)2 (5.9723522 x 10°!) = 


and k = 107.0731376 and k* = 11),6.6568 
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TABLE 1 
COMPUTATION OF VOLUME, DENSITY, AND MASS OF THE EARTHs COMPUTATION OF k 


DIMENSIONS OF THE INTERNATIONAL ELLIPSOID 


CX = 1/297 1+ ¢ 2 = 1.0067681702 = 

a = a meters 23 0005568098” ee 
b = wae. meters poe Oe 

e. = 0.081991889 we Se = = gs 

e? = 0,0067226700 ere nae tina 

e , = 0.082268889 == 

5 = 0.0067681702 } i i aid 

e°~ = 0.0005568098 sop ere = 020817158230 


(<= (2. x 3. 15926595)" 6,378,38800. = 3.3805);56),00 


V = h/3 TT a2 = 1.083320781 x 1027 cm 





ame 358 mass -— |) 


M= VO n= 5.972352200 x 10°! gm 


Pe a a 3 : (86,00)?(5,972352 x 10°! 
ia em cm? gn. 
6 gece G eg a Gs (6,378, 38800)? 


107 20731376 





. DETERMINATION OF THE ORBIT 

Since the arbit of the satellite is defined by six parameters then 
six independent quantities mst be given by the observations in order 
that the parameters may be determined. A single observation gives two 
Quantities, the angular coordinates of the body. Therefore, three 
complete observations are sufficient to define the orbit. [byP.219 

The apparent positions of the observed satellite are obtained 
by measuring its angular distances and directions from a background 
of stars on the film of the Baker-Nunn camera. The stars on the film 
can be identified from catalogues containing their right ascension and 
declination. With these coordinates known, then the right ascension 
and declination of the satellite cante computed. Catalogues containing 
the observed values of time, right ascension, and declination reduced 
from the camera film are published about six months after the observa= 
tions were made. These reduced values are quite accurate. However, a 
more precise reduction of the films is made by the Smithsonian Institu- 
tion Astrophysical Observatory Photoreduction Division. These results 
are published and made available about 18 months after the observations 
are takene 

In this paper the computation of the orbit of Explorer I (1958 Alpha) 
satellite will be reduced from observations of right ascension and 
declination from the Baker-Nunn Camera Tracking Station at Olifants= 
fontein, South Africa. Date of the observations was 18 December 1960. 


The observed data [7, p10] 


Time (UT) ) oe 6 








225 og” 308.24 gh o7™ 65 — 7° 22'00" 
22h 317 obS.h2 iy 28" SiS — 52° 06'00" 
22h 321 938.10 8" 337 oo® —}9° 25'00 


These data are not those precisely reduced by the Photo-reduction 
20 
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Division but the precision of the observations is as follows: 4 Peo 


Standard Error in timing ( J+ ) co5S= = =Ge = woe? 
Standard Error in direction ( 9G, ) 2'.7 < 0, = 


The geographical coordinates and height of the tracking station are 
[8, p.7]: 


D 
H 


There are two classical methods available for determining the orbit— 


25-257! 3s 70S 
28° 1h' 51",.10 EF 
15), meters 


A 


i 


one developed by Laplace in 1780 and the other by Gauss in 1801. Both 
methods were devised for use in determining the orbits of comets and 
planets about the sun but are also applicable for use with artificial 
earth satellites. Since the development of these two methods many 
papers have been written on the theory of orbital determination but 
very little that is theoretically important or radically different 
has been added to the methods of these two men. The Laplacian method 
is not as complicated as the one devised by Gauss but it has two weak- 
nesses. It is inaccurate when the intervals between observations are 
short and it assumes that the observer is in orbit about the focus of 
the ellipse. The Gaussian method, even though it is a approximate 
method was chosen to determine the orbit. There are other more modern 
solutions available but these are not as accurate. 

The first step in the calculations is the conversion of the 
Universal time of the observations to sidereal time. This is required 
in order to introduce the vernal equinox into the coordinate system. 
The sidereal time at any instant is the hour angle of the vernal 
equinox. 

The origin of the coordinate system is at the center of mass of 


the ellipsoid. The x -axis points toward the vernal equinox, the y-axis 





[=a = -« 
— _. = 


a ¢ <= 





| ‘t= = 


~~ — = « 


| > — en 


— _ = 


=f. = «=. ©< = 
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is directed toward a point 90° eastward from the vernal equinox, and 
the z-axis coincides withthe north pole of rotation. The x-y plane 
lies in the plane of the equator. 
The rectangular space coordinates of the observer are computed 
fromthe given latitude nd longitude of the tracking station. Then, 
by changing the signs of these values the coordingtes of ue: conve of 


the earth are obtained. The observer's position is computed as follows: 


x43 = (N+H) cos D cos Oi, (3 ="dy (2; 3) 
yj = (N+H) cos D sin 0 1; 
2j = (N+H) sin@ (1 - e?) 


The position of the observer, the earth's center, and eventually the 
position of the satellite are all converted to rectangular space 
coordinates to provide a more convenient and a common system of 
coordinates. 

The next step is to compute the direction cosines from the 
observer's position. These are functions of the observed (apparent) 
right ascension and declination and are obtained from: 


cos o, cos KX , (j =81, 2503) 


f. , | 
= b : 

84 = cos J sin 

hj = sin 6 j 
The orientation of these axes: f=<axis points toward the vernal equinox, 
g-axis perpendicular to f<axis, and the h-axis parallel to the axis of 
rotation of the earth and perpendicular to the f~ and g-axes. The direc- 
tion cosines give us a line in space fromthe observer's position to the 


satellite. The computations of these three steps for the three observa~ 


tions are contained in Table 2. 


ae <— = « 
~—™ + -_ 


= oe 
a, ° — 





TABLE- 2 23 


COMPUTATION OF SIDEREAL TIME, COORDINATES OF EARTH CENTER, DIRECTION COSINES 


Date 18 December 1960 18 December 1960 18 Dec. 1960 
Sidereal Time At O° UT Shp bu25 1:29 Cel 6-25 04419 56-25 6419 
UT Time Interval 22=29-10.140 22=31—0)) 1420 22~32-03 .100 
UT To Sidereal Time 3=)1.633 J—1 9116 3~),2 «107 
Greenwich Sidereal Time 20—22—10.020 


~L7—™ a 2 20—P1=Lle BS 
(Add East, Subtract West) +1-52-59.1,07 + 1-52-59 1,07 +1-52-59 0,07 


Local Sidereal Time (19 Dec) “Go TTS192 303 
g 930)! 08" .985 93932'h7".880 93,7130" 495 

sin 04 9985656) 998081178 .9978109), 
cos Oy = 05354139 ~.06186091 - 06613118 
ao sh_o7m_o68 h_ogm_c),8 gh_33_908 
Ct’ 76°h6! 30" 112013130" 128°15t00" 

sinCX 097347917 092570563 078531693 

cos(k 022877565 0378247 = 61909395 
& 4722100" § 516100" § 49925100" g 
aes = 73570317 = 677824315 =6759),6058 
cos () 67730410 262796306 ~65055333 

ee 25°5713".7 S oe 

_— + ~ 013773790 <__-— 

By —————— 01916117 <———_—— 

io ——: 089910263 eee 

a 6,378,388.0000  <——— 

e@ ————S 000672267 <——__—_ 

L- ee —— »99327733 = 

e*sin’D ee .00128816 Se 

1 = e*sin2® , ee 09987118 <a 
(lee@sin?Q )2 ae 99935571 a 


a/(1-e%sin® PD \z Nv o--— 6,382,500.18 meters <—_— 


es 154 meters <——— 

sie —— ~0002),207 ae 

N/a —__—»> 1.0006h171 oe 

N/a + r = N+H ea. 1.00088678 oe ai 
-(N+H) cos D - mG =x +.0),818189 + .05566863 + .059511hh 
-(N+H)cos =y ~-.89860916 -.898176):3 - «89793000 
~(N+H)sin (1-02) =z +.13518070 + .1;3518070 + .:3518070 
- 2 awe cu + .181)95069 ~ 23752372 ~ 40275363 
cos() sin(y! =g¢ + 66593143 +.5813089) +.5108905) 
sin 6) =h = 673570317 07782315 ~ 075946058 








2h 
The next step is to determine the distance of the satellite from the 
observer's position and from the center of the earth. This is done by 
solving the triangle formed at the time of the middle observation between 
the center of the earth, the observer, and the satellite. To solve the 


triangle the satellite is "stopped" in its orbit at the time of the 


second or middle observation, (Tp). 


qy 





Figure h. 
The geometric situation is depicted in Figure h. Since SOE form a 


triangle at To the quamtities ro and Ro must satisfy the equation 


luyp-201] 





RS =r? + oe - er (0) cos \/ 


This equation treats the geometric situation. In the Gaussian method 


a formula relating the two wnknown quantities r and & applies dynamically 


[4sp-235| é 2 
as 
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where kK, A » and Kj are determinants involving the coordinates of the 
earth's center and the direction cosines from the observer's position. 
These determinants and their values are: 

ji <; > 3 2 = — 0.0020)368 
hj h2. h3 
ff, =X] + X3 - 2x2 £3 


1 
K = -- | g y,4+ 93 - 2y2 
ha zt + a ~ 229 h3 


+ 0.00018897 


0a 
WJ 
8 


f X1 + X3 f3 a 
Ki= - gl J+ ¥3 g3| = 045716432 
h, 21+ 23 h3 


The solutions of these determinants are shown in Table 3. 
The two unknownsry and Ro are expressible in terms of a single 
unknown, the angle @ » by means of the equations [hp -235] : 
P. cos Y 5 = x0 fp + Yo ep + 22 ky (0 = Yaa 
P> sin’, = + .48507360 
Po cosy =-K_ = ~.781551K6 


The computed values for Nhss N, and m are: 


N sin m 


Nicos m 


Yo = 150° S8t 12",06) 
N= * 9.9198)731 
m = 148° lot 26",395 


The solutions for these quantities are contained in Table he 
N, m, and M are intermediate parameters used in the solution for the 
angle p » where ls p.212| 3 
hi A a N si 3 
rs sae = + 0.163)1)533 
$2 Kl] ) 
The sign of N is chosen so that M shall be positive. 
The equation for the solution of the angle D is bh p.213| : 
sin + = Msin (@ + m) 
which is a transcendental equation. The quantities M and m are know and 


Mis positive. One solution of the problem results in the position of 














<—“itené . - - ‘ 
— aoe = = 
i Cy Sm, 
. =—@ 5 





° Sein oo — > =" 
a 














_ 


Ky 


a (yl + a 


(at X3 = 


2) 


TABIE 3 


SOLUTION OF DETERMINANTS 


ry 2 83 2 
£3 g1 he 
£3 89 hy 
fi g3 ho 
fo g1 3 


Yeats 


£36 (24+ 23-222 
Y1.+ ¥3-2y2) hy 
1 g3(21+ 23-222) 
(x1 + X3-2x2) g)h3 


fy (y1 +3) 


h3 


(x1+x3) 83 hy 
£3 g1 (21 +23) 
£3 ae y3) I 


(2 83 eit Z 
(x34 X3) 21 


2 


~ .068 107 82 


0089 276 57 
206 66) 15 
0172 245 98 


0118 938 Sh 


e000 O21 92 
001 369 62 
0 


-.000 055 20 
0 
e001 82), 68 


2211 hil 83 
—00 478 02 
2068 900 28 
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TABIE ), 
SOLUTION or, N, m 


e 


x9 fo =  =4@13 222 62 
zo ho s -.338 676 ho 
X2 f5 + yO Bo+ 42 ho = = O07 O17 01 = P, cos Yo 
= 2003 099 00 
= — 720 ah 
z = 0189 382 2 | 
xh + yo + 35, = 3 2 
(x§+y + aBys = ~999 600 99 = P, 
xof>o + ¥280 + Zoho/Po2= - 287) 365 89 = cos Y 2 
bias = 150° 58112" 06) 
sin YS = 485 267 23 
Po sin Ns = «485 073 60 = N sin M 
P2cos Yo--K = -.781 551 h6 = NcosM 
N¢ sin? m = «235 296 0 
N@ cos2 m = .610 822 68 
N2(sin@m+ cos2m) = ‘ 
N = 0919 8h7 31 
Pp. sin Vo /s = 0527 3hl 43 = sin 
ia cos Wak A = -.8),9 653 h7 = cos m 
m = 148° 10'26".395 
A = -,002 03 68 
P> =  .998 803 h5 
sin? Y = ell 272 81 
N = 0919 847 31 
ts = 011-485 88 
Kl 3 = —57 16h 32 
hA P3NstwY = -.000 858 2h 
t® Ky =  =.005 250 93 
bh A 03 Ning = 163 Ws 33 = i 
te KL 


k (T - T1) 


Solution for t: t 


T3 = 22h 32M 035.10 
T = 22h oo 108.1) 
T = Eat = e 
- T, = .002 001 85 days 
k = 107. 073 1376 
k(T3 - T1) = 0.21) 3h) 36 = 2t 
4 = @.107 172 16 
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the observer and is not applicable. That is, ifr=0, J=R. The moper 
solution for the angle PD must satisfy the inequality lh, p.213| : 

OP-1-¥ 
In our problen, Y equals 150° 58'12".06) so p must be less than 


2°02'. The solution may be obtained graphically be plotting the curves 


= sint 


yg = M sin (DP +m) 


Once the approximate value is found from the intersection of the two 


|b, p.213| $ 


curves then the exact solution may be found by numerical trials or by a 
differential form la. The solution for the angle is lengthy and time - 
consuming. The approximate value from the intersection of the curves on 
the graph was found tote approximately 23°. The exact value, solved for 
by numerical trials was: 
D = 23° 15" hon,16 

With this value known we can now solve for the two unknown sides of the 
triangle SOE in Figure h, which will yield the distance of the satellite 
from the observer's position and the distance of the satellite from the 
earth's center for the second observation. The formulae used and the 


results are [a De 255| 


Rp = Pp» eae 2 = 1.228 150 03 


sin 


Po = P, ng y+) )= 0.25, 280 87 
sin 


The observer - satellite distances for the first and third 


observations are given by Ly De 255] : 


fi f3 - (x1f3| _(1,3)[xof yah xaf ee fof 
Bel» ae 3 (3,3) 3283) 3 3a + age.) 





ay 


f)f3 ey f)x7/_ (153 £ix2 [4 £4x3 ot fae 
Be | *3 a 9¢/ |B1Ye2 "| ely; = 9 | sue] 


where r 
$434 = a = 2.93 11079 
ae sides = 1.943 20695 
L+ Et oR 
+o 4 
(3) = aie = 045161323 
= + a 
G3) z — - 1.50993233 
’ dats $< - 
where 
42 KU =  0.10717218 


et Ta) = 0.035185 


q 


The time of the second observation, To, is the origin of time. Upon 
substituting these values into the above formulae and solving the 
determinants, the observer - satellite distances for the first and 
third observations are found to be: 
ry = 0.23531756 
r3 = 0.27239982 
The detailed solutions for Ro, rj, r2, r3 and q are listed in Table 5. 
The next step is to determine the geocentric space coordinates of 
the satellite for the time of each observation as follows lh p.256| : 
r 0223531756 0.25),28087 0227239982 
X=rf-x ~0.01171927 -0.11606637 -0.16922145 
YS=rg-y + 1.05376378 +1.04599217  +1.037096h9 





TABLE 5 
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DETERMINATION OF Ro, rg, 4s Fly 3 


Po sin Be, = 04185 073 60 
= 0.39 962 82 
ps4 ss Ne an Ss 1.228 15003 = Ro 
W = 150° 58112" .06 
a xz 0 1,9" 16 
Sa aia? 1OL" 22 
feo D ) =  0,10047158 
se sin CWE + ® ) = 0.100319 
Po sin( (Yo + @ )sin = 25128087 = 1ro 
T3 32™ 035.10 Tg = 22> 31° O4S.2 
To » 3s 228 31" oh .2 T] = 22% 29M 108.1) 
%-T = ee tee: 585.68 T> -T; = (-) 1” 548.28 
To = T] = (=) 1 55.28 
T3 = To = (+) 585 .68 
T] + % = (-) 558.60 
T] + T3 (days) = -0.0006)352 
k = 107.0731376 
k (T] +73) = =0.06890371 
k (T]+13)/2 = -0.03445185 =i 
t  0.10717218 Ro? 764099326) 
26 = 0. 21,3 4,36 q -0.0345185 
t2  =0.011,8588 q/2t -0.16073131 
R 1.22815003 a/t., -0.321)6262 
R2 1.852),8316 +2 /2Ro: 0.00310013 
Qo? —-3. 7096632 t2/lRo> —-0.00155007 
1-a/t tel = 1.32456275  1+q/t + +2/aR2? = 0.68163751 
ro, 2, = 4 bon : = 2.93h11079 
2+ mE + =? 3, 
t vs = (2) = 1.913206 
1-8 Si /1+ 4s oe m } o9h3 95 
2 
q. +2 /i et, = [lade 09,6111)6132 
adage / +A ae * RB = osuswe 
2 
/ b-A + y= (2) = 1.50993223 
a > 








a 
_ 6s os) 

7.) 

- 

-_—-F"¥ eo 


TABLE 5 (Continued) 


X71 e 
Lae 
x183 - yif3 =| | = 
X983 
yet3 _ | x@els} _ 

X983 = Yof3 ~ | Y23) ~ i 
23) xof3 - > 
\és Y283 
X38 
y3f3 > x3f3 = 

=3ep8 3 | RS ‘ 
a pees sy <2 

9 9383 
3283 = 
52*3 Sigg) 2 
fo83- Bots = ee = 
1 fof3 | < 

Be 9 ea 
“183 
8] 3 fl f a - on 

#183" 8143 * 4 : 
£1 f3 
6183 aa = 
ry = = .08111749/.34471L99 
fiyz te 
81) 

fiyi- ax, = sap | = 
7 £41 = a 

9¢ B1Y1 
£19 , = 
g1x9 fix2; _ 

tyyo ~-enX> = ae = ae 


L ) 2 
335 E192 


0.02)61567 
0.36191810 


@ 


0.028),),058 
063617382 


02333350324 


0.97794863 


0.030),0383 
0 3616),);57 


Oo33Len074 


0.64,366931 


0.1213),862 
0.231229 


0011277507 


0.08)11,060 


0.0791628), 
0226555215 
0.3uh71H99 


0.081117h9 
= 0.23531756 


0.13924011 
0.03176832 


0.17 100043 


0.08800320 


0.13917306 
020367063 


0017587769 


0.26556339 
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TABLE 5 (Continued) 


f1y3 - 0.13913)87 
81x3 _ f£yx 7 0.03923836 
fiy3 - 6x3 «= aa - - 0.17837323 
£1 £0 0.09007)22 
gifs _ [ £y£5 - 0.15660923 
82.5 81f2 = | oes] = = 0.21,668315 
a2 ay cae = 0.0947133k 
cs - = 0409390030 
8183 3 
=  .09390030/ 63447199 = 027239982 


tis = 
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7 =rh-z ~ 0.60830457 ~ 0.63307305 ~~ 0.64205763 
R =\(x2+y*+ 2°) 1.21679h93 1.22815003 1.231139}0 

There is a check for accuracy of the computations up to this point 
in this last step. Note the agreement in the value for Ro computed from 
the space coordinates and that determined from the formla, {/> sin Y 9)/ 
sin PD » in Table 5. 

The inclination of the orbital plane and the intersection of the 
orbital plane with the plane of the equator (XY plane) can now be solved 
using the geocentric space coordinates in a ratio of determinants. This 


solution will produce the first two parameters of the orbit that we are 
seeking: i and (2 lus p.256| : 

2419) ; YO 

Loko Xo Yo 


e e = Yy 27 
See Table 6 for the evaluation of the ratio of these determinants and 


the detailed solutions for:-i and () e 


Then, 
cosi= <Q i = 32° 31,820" 25 
ae + Be + C2 
sin ( sin i= A 


>) /a2 + Be + c? 


ge sin i = =B 


(2 = 206° 00127",603 
ake ieee ce 


If the inclination was greater than 90° the sign of the previous 
two equations would be reversed. The values i and ¢. uniquely determine 
the position of the orbital plane in relation to the equatorial plane 
of the earth [h, pelh6|. 

The next step is to find the auxiliary parameter u for each observa= 


tion. Solution of this value is necessary in order to find the "longitude" 





TABLE 6 


3 


EVALUATION OF DETERMINANTS A: Bs C 3 SOLUTION For a anp $2. 


YZ 9 
YoZ 
2°1 Y4Z 

Ti¢2t= elie” ae 

ZyXo 

Zox ZX} 
Z 1X9 - ZoXy = Loko 

X1Yo 


Xp¥ 
2X} XY 
X1¥o°- Xo = as | 


> /a2 + B2+c@ 
2 ee Ee 
C/\/ac+Bo + C 


i 
sin i 


4/- [x2 + Be + 6? 
ee oe 


Ao 4 Bo + #/sin i 


of [a%4B2 + 62 


B 


- /a2 + 3? + c2 sin 


(2 


0 .667109))5 
-0.63628182 
—0.03002703 


0.07060370 
0.0071915 
0,003 
=0 01225826 


=—Q ° 12 23 065), 
01100828 


0213058811 


084271286 


3293120" .2h5 


0253836329 
—0..23606766 


= 043 84,9138 


0 4838612 


=) 89873535 


206°00! 27" .603 


cos il 


cos $2 
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of the satellite from the position of nodal passage. The longitude of 


the satellite from the node is called the argument of the latitude. 





Figure 5. 


Refer to Figure 5. The formlae that define uj, uo, u3 are [uyp.236] : 
cos us = cos bs cos (a, - 82 ) (j = 1, 2, 3) 
sin uj cosy = cos bj sin (1,5 - (? ) 
sin Us = sin by / sins 
To determine the parameters a, e€, and w,we mst solve for the ratios 
of the triangles formed at the time of each observation. With the geo- 
centric coordinates and the range of the satellite from the earth's 
center known the above equations may be substituted with the following 
expressions ld, p.253,257| ° 


eet Sa j = (1, 2, 3) 


cos i sin uj 


2 if oe? = * =? 


Bt 


a ear? + + can. 


cos us 


Then for abbreviations 


OC, 233 = ug = 2°obh197583 


ao = uz - uz = 7°45516083 


0.04262050 radians 


0.13011710 radians 


O, = up = uy = 5°.01318500 = 0.087)9658 radians 
Then the ellipse parameter p is determined and used in place of the semi- 


major axis, ao 
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The equation that defines p is lh, p.257 | 2 
kVp! (T3-T]) = RS Op} + R,? On(2 S 3-Gy)+ RS Oo(2 0,- 03) 
601 03 6 03 6 Oy 
Using this formula the value for p was found to be; 
p= 0.8305535h 
With the values of p and u now known the eccentricity, e, and the 


argument of the perigee, w, can be found from [us p.257| 


e sin (ul -w) = R3(p-R}) cos J» — Ry (p-k 3) 


e cos (u,;-w) =p — Rj 
L 


Then the semi-major axis of the orbital ellipse is defined by 
2: a. 
The computed values for these three orbital elements are: 
e = 0.32815828 
a= 0.93078812 
w= 82°5)')2" 337 

The detailed computation of these parameters and the auxiliary para=- 
meters uj, U2, u3, and p are listed in Tables 7, 8, and 9. 

The next and final orbital parameter to be computed is the time of 
epoch, To (time of perigee passage). This value is computed from Kepler's 
equation: 

M=EHE-esink 
after the true and eccentric anomalies are determined for the three 
observations. Values for To are? 


Observation #1 Observation #2 Observation #3 
December 18.915)9)2 18 .91,53182 18 .91455056 


Detailed computations of Tp are listed in Table 10. 


To verify the results of the computations, the orbital parameters 


were determined again using data published in Smithsonian Institution 
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Astrophysical Observatory Special Report #78 dated 25 October 1961. This 
publication contains both the mean elements and the smooth elements of 
Explorer I (1958 Alpha) satellite based on 82 observations made from the 
various Baker-Nunn Tracking Stations for the period 16 Decémber through 
29 December 1960. The published smooth elements are: 

To = 37291.0 Modified Julian Days (Smithsonian Days) 

(10°.6 +2)+7°.275 (t-To) + 0°.000427(t-T,)*+ 0°.22h cos w 


W 


0 


(1869.278 + 9)=°.896(t-T,)-0°.000298 (t=Fo)* + 0°.00887 cos w 
(33°.210 +3)=0°.00055))(t-T)+ 0°.50x1079(t-1,)2 -0°.00h73 sin w 


. 
ai 


= (.09971 +8) ~0.000202(t-T,) + 0.30 x10~7 (t=1,)°+ 0.00045519 sin w 


M = (.9016 +8) (13.43513h +5) (t-To)+(0.0002839 +6) (t-To)@ 
=(.11 +1) x1075(t-T,)3 -(.000631 x1073) cos w 


© 
t 


Standard error of one observation: m= 1451,)3. 
The angles are expressed in degrees, e in terms of the earth's semi-major 
axis (6,378,388 m.), and M in revolutions. The Smithsonian Day for 
December 0, 1960 is 37268.0. Using the time of epoch, To = December 
18.9153182 from the middle observation of our calculations and redesigna= 


ting it as t, the term (t = T,) will have the following value: 


Date Smithsonian Day 
December 0, 1960 37268 .0 
Date of Observation +18.911,53182 
Date of Observation (t) 37206.91 53102 
Epoch — (To) 37291.0 
(t -To) =)1.00 Ss OOLO days 


Detailed computations of the elements are contained in Table ll. The 


computed and published parameters are: 


Published Computed 
i 33°13100.)02 3203)'20.2h5 
0) 206°16' 3" .067 206°00!27" .603 
e 0.10011028 0.32815828 
a 1217098895 0.93078812 
w 120°57 #53" 675) 82°5)! 2" .337 
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COMPUTATION OF THE AUXILIARY PARAMETER, u 
Observation 1 Observation 2 Observation 3 


Z 372 —0 049992366 ~0.515)6882  -=0.52138792 
sin 1 053836329 0.53836329 0453836329 
(Z;/Rj)/sin i = sin u -0.92859909 -0.9574739  -0.96846856 
is 0286601592 0285168110 0.8)218232 
1/85 ~0 00963126 -0.09)5050)  -0.137)1760 
s ——> -0.1,38),9138 << — 
cos ——— ~0.89873535 <————— 
Y/R sin ~0.3797052 -0.37345482  -0.36928969 
X/R cos 0.00865595 0.08h93502 0.12350205 
y/Rsin ()+ X/R cos (2 = cos u -0.37108K57  -0.28851980 -0.2457876) 
vu 218913102" .845 253913'50".311 255°)0821" oh 
gun = Oy 2026 31" 113 
oa (radigns) 0.0262050 
u3z-U, = 7927918" 579 
Jo (radians) 0130121710 
Upg=U) = J3 5°001 7" 66 


Gs Gadians) 0.087),9658 





TABLE 8 a7 


DETERMINATION OF AUXILIARY PARAMETER, p.- 


kVp (13-2) = oe 
601 a 


BR 8 
Ros 1 O36 %, oF 


& gs? Oy) 

Rae gat O3- F1)/6 C3 

Gi = O03) 

ef de 4 O1- Sa Oy 
art) 

' ee = 


ape whet) zaVp = 


+%+e% = 


Ri° P22 oa On): R3° P22 a J) 


Oy 


1.148058990 
1.50835250 
105164300 
0,0022029); 
0,00332281 
0.02237h90 
0.185061 
0219265006 
0.13237266 
0.529798 
0.0)857637 * 
0.19731517 

=( »00225558 
0225572300 

=0.0017h0h0 + 
051953211 
06214336 
0.9113h710 
0.83055354 





TABLE 9 


COMPUTATION OF THE ORBITAL PARAMETERS €y a, We 


e sin (uj-v) = R3 (p-R1)cos Oo-Ry (p-R3) 
RZ ' R3 sin ) 


10231h39)0 
es 0 .3862139 
sy: > 0.9915),671 
1.21679h,93 
Rs ) =0.1,0088586 
R3 (p-R, ) =0.147563287 
R3(p-R1)cos oF =0 17161221 
Ry (p=R3 ) -0.18779588 
Ry (p-Rj )eos Oi o-Ri(p-R3) 0.01618367 
sin Oo 0012975025 
RUR3 sin OF 0.19441897 
e sin (vu) - f) 0.08321,121 

e cos (ul = w) = (p-R1)/R1 = =0.317)2521 
ensin, (uy) 000692910 
e°cos °(ul-) 0.10075876 
- |sin (uj-" }4+c0s?(uy-W)| 010768786 
0.32815828 
ns = a 0.93078812 
e sin (u 17) 0,.0832)121 
sin (u ,-7) 0.25366177 
e cos (uj-w) 0 ..317),2520 
cos (uj-w) =0.967 29298 

(uw) 165018' 20" 508 
U1 21:8913°02" .Bh5 


w 8205),'h2" .337 














V = u-w 
v/2 


l=<e 

l+e 
(l-e)/(l+e) 
Vil-e)/(ite 
tan V/2 

tan E/2 


E/2 
E 


E (radians) 
sin E 

e sin E 

M (radians) 
M 


k 
ae 


Vad 


n 
E-e sin E/n=xt-T, 


t3 (days) 


TABIE 10 
COMPUTATION OF TIME OF EPOCH, To 


Observation #1 


165°18! 20" ,508 
83939910" ,.25) 


Observation #2 


170°19!07" 97h 
85°09#33" ,987 


tan B/2 = V(i-e)/(1+e) tan V/2 


——__—> 0.6718)172 
—_—_—__> 1.32815828 
———— 0.5058)),62 
> 0.7112275) 
707555753 11.808});28 
5 25159787h 8 .398)7920 


79°93 #27" 723 
159926155" .)))6 


83°12#35".307 
166025?10" .61) 


M=- E-e sinE 


Observation #3 


17295139" 087 
86922 #),9" 5443 


———L— 
x 
ak 
——— 


15.80837), 
11..2)335095 
84 55#02" 65), 
169°50!05" .308 


2 .96),17658 
0.176);8680 
005791560 
2.90626098 
166°930'59" , 358 


—————————— 
—— 


SS 


1192351450 rad./day <—_—-—— 


27829051 2 69057022 
0.3510)4533 0.23),80935 
0.125198h3 0.07705)63 

2 .66770698 2 082751559 
153925?16".711 162900!16" .966 

n = k/yVya3 

— 107 60731376 
naan 0.806)0366 
——<$__> 0.89799981 
0.02237350 002371378 


j = (1,2,3) Dec 18.93692292 


To 


Dec 18.9145)9h2 


Dec 18.9382);560 
Dec 18.9153182 


0.023720 


Dec 18.93892),768 
Dec 18.91:55056 


TABIE 11 2 


ORBITAL ELEMENTS DETERMINED FROM PUBLISHED DATA 


(t=T,) = ~); .085),6818 days ws 110° 8 § ho " o 896 
sin w= 0.9375525 
(t<To)** 16.69105025 cos w= —0.35529231 


(t-T 9 )3=— —68.19075)69 

wo = (1409.6) + (7° 0275) (ly 008546818) + (9.000427 ) (1669105025) 
@(° 22h) (=.35529231) 

= 110° 57! 53".75 


= (186°.278)—(° 0896) (=; 08546818 )—(°.000298) (1669105025) 
=(°.00887 ) (=.35529281) 


= a 


206° 16! )3".067 


i = (33°.210)=(°.00055)) (=l} 008516818 ) + (°.000005 ) (16269105025) 
+ (9.0073) (93475525) 


= 33° 13" 00.102 


= (.09971)=( .000202 ) (=l.085)6818 ) + ( .00000003 ) (16.69105025 ) 
=( .00045519) (093475525 ) 


fj. 


o 


e = 0.10011028 
a = q/l-e = 1.053706/.89988972 
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Comparing the above results it is apparent that the position of the 
orbital plane is correct but that there is a discrepancy in the orienta- 
tion of the orbit in its plame and in the size and shape of the ellipti- 
cal orbit. Reviewing the computations, it can be seen that the auxiliary 
parameter p, which has a value of approximately 0.830, is too small. In 
the case of an elliptical orbit it should be slightly greater than one. 
Since p is a function of the auxiliary parameters uj, it is believed that 
the source of the discrepancy rests here. This would affect the results 
in the computations for G; and for the argument of the perigee, w, 
also. To avoid continuing the problem with this discrepancy present I 
decided to use the published orbital data from the 82 observations rather 
than to continue with the non-realistic values for e and a derived from 
the computations. The orbital parameters computed for the times of 
observations in the next chapter are based on those published in Smith- 


sonian Astrophysical Observatory Special Report #78. 





Se DETERMINATION OF THE GEOGRAPHIC COORDINATES 
OF THE UNKNOWN STATION. 

There are several methods available to determine the position of an 
unknown station from observations of a celestial body. Two of these 
methods were devised by William Markowitz of the U. S. Naval Observatory 
for his well known Moon Camera Method and can be adopted for use with 
artificial earth satellites. These methods are called the linear, 
which gives the coordinates of an unknowm station, and the differential, 
which provides a way to find the corrections to an assumed approximate 
position. The differential method has been selected for use in this 
particular problem. As mentioned earlier in this paper, only two 
approximations in the differential method will be made. 

The procedure to determine the coordinates of the unknown station 
is as follows: 

1. The orbital elements of the satellite are known for a particular 
time, To, the time of epoch. Now assuming these elements are COneEt 
the orbital elements are then determined for the times of the observa~ 
tions from the unknown station. 

2e Using these elements the geocentric right ascension, geocen=~ 
tric declination, and the range of the satellite from the earth's center, 
will be predicted for the times of the observations from the unknown 
station. 

3. The differences in the computed coordinates and the observed 
coordinates of the satellite will furnish us with a set of coefficients 
for three differential equations tobe solved for each observation. i 
solution of the differential equations will produce the increments A CL 


t 
and A 6 which are added to the approximate coordinates. Using these 
bly 
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new coordinates, the differential process will be repeated and the resul= 
tant increments will be added to the approximate coordinates again. This 
will be the final result. 

Refer to Figure 6. 

The observer, 0, views the satellite along the line OS? at a distance 
of re The coordinates observed from this position are the apparent right 
ascension, x » and the apparent declination, 6 » From the center of 
the earth, E, the satellite would be seen along the line ES! at a range 
of R. The coordinates observed from this position are the geocentric 
right ascension, CX , and the geocentric declination, 6 o The coordi- 
nates of the satellite derived from the orbital elements are the geo= 
centric coordinates and will differ from the observed coordinates determined 
from the camera film by the amount of arc, SS', in the orbital plane. 
This difference enables us to determine the observer's position on the 
surface of the earth. This difference between the computed and the 
observed positions of the satellite should be very small but not so 
small that the two positions would be at the zenith of the observer. At 
this position both the computed and observed positions would have the 
same value and the observer's position could not be determined. Another 
indeterminate situation would exist if the satellite was observed at 
the horizon. Here the parallactic angle would be too large and the 
effect of refraction would make the observations useless. lae, Ped ° 

Proceeding with the problem, the following data are known: 

Location of Observer: Curacao, NWI. 

Approximate Position: p # 12° 13'23".21 N 
68° 55:07" ,.86 W 

0235199282 
0.913105082 


va. 
Za = 0.210569117 
Satellite Observed: Explorer I (1958 Alpha) 


6 8 ft 


Xa 
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Universal Time of Observations: 23%5)™,7S.o0), 


23h51Ms5),8 .98 
a3heom 05 61,5 


The approximate rectangular space coordinates, Xa, Ya» 2a, are found in 
the same manner as that used to find the coordinates of the known station. 
This procedure is discussed in Chapter | and will not be shom here. 

Using the above data, the orbital parameters for the times of obser= 
vations are computed. These computations are contained in Table 12, 

With the orbital parameters know, the next step is to compute 
the eccentric and true anomalies of the satellite using this formla 
(usp-262] for the eccentric anomaly: 

E = M+e sin My_S° sin 2 M+ 


oo080 0606 


where Mis in radians. The additional terms in this formula are not 
listed or used in the determination of E. Moulton states that because 
of the small eccentricity of the orbit the series converges very 
rapidly and the first three terms give an accurate value of . 


Then the true anomaly V is computed from the formula: 


tan bv = Viite)/(1-0) tan & 


The computations of the anomalies for the three observations are listed 
in Table 13. 

The next step is to determine the geocentric right ascension, geo= 
centric declination, and geocentric distance of the satellite for each 
observation using the formulae: 

sin O = sin i sin (w+y) 
MX = (? + arctan [ cos i tan (w + v)| 
R =a (1 =e cos E) 


where a is the semi-major axis of the elliptical orbit in terms of the 
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OBSERVATION 1 OBSERVATION 2 OBSERVATION 3 


Location uhes ao,N oI. Curacao oN WoT Curacao yN.W.To 
Time of Observation (U.T.) 3 175 ol 23%51%5),5 .98 23%50m7 95 1150 
Date of Observation 18 Decanber 160 18 December 1%0 18 Decanter 1960 
Smithsonian Date of Observation(t) 37286.99h29hli) 37286.99)38G)3 37286.99)586l2 
Smithsonian Date of Epoch (T,) 3729120 37291.0 37291.0 
(t-T,) 5005705556 ~~ 00501357 aloe 
w (Mean) 111923827" ,7255 112923130",1299 = 111923'35" 12 00 
sin w 093111290 © 93110865 093109929 
cos W = 6 36473110 = 03674195 - .36),76583 
(t=To)* 16.0)567700 16 .0)491077 16 604332095 
(t-To)? 64027425751. ~-6);.26983387  =61,.26010156 
w= (1409.6 +2)+7°.275(t-To)+ 0°..000h27 (t-To + 0°. 22h) cos W 
7° .275(t=To) ~299 01150792  ~299.1083935  -29°.13936836 
0° .000427(t-T,)2 0° ,00685150 0°,00685119 0° .00685050 
0°.22h cos w 0° 08169977 =00,08170220 ~0° 08170755 
Wo 140°.6 +2 a” 66) eee 140°.6 + 2 
Wy 111°.3836381 § 111°.383096) 111°.38577h59 
Wi, 1119 23901".118 122023t03'".515 121°23'08", 789 
(ds (1869.278 +9) - (°:896(t=To)=-0°.000298(t-T,.)“+ 0°.00887 cos w 

4°.896(t-T,) 19°,611934h0  19°.611),8h46 19°.610)9)50 
0°.000298(t-T,)2 ~00,0078161  =-0°,00)78139 ~0°.00h78091 
0°,00887 cos w ~0° .00323516 —0°,00323526  -0°.003235h7 

To 186°.278 +9  186°.278 +9 186°.278 +9 

t 2050.88191763  205°.881)6781 205°.8807812 

t 205°528Sh" 90B =. 205952153. §=20595211,9"721 


i = (339.210 +3) -0°.0005Sh(t-T,)+ 0°. 50x1L0-?(t<I'9)2~0°.00473 sin w 


0°.00055(t-To) 0°,.00221916 0°,00221911  0°.00221900 
0° .000005(t-T,)@ 0°,00008023 0°.,00008022 0°.00008022 
0°.00473 sin w ~09, 0044016 -0°.00hLOL) =0° ,00)0)10 
i To 33°.210 + 3 339.210 7 3 339,210 +3 
it 330.20789523  33°.20789519  33°.20789512 
it 33°12! 281,1:28 33°12128")23 33°912'28"1e0 
e = (.09971+ 8)-0.000202(t=To) + 0630x1077 (t-Tp)*+4+ 0.0005519 sin w 
0.000202(t-To) 0.00080915 0.00080913 000080909 
0.00000003(t-T,)* 0,.000000),8 0.000000),8 0.000000),8 
0.000),5519 sin w 0.0002383 0.000);2383 0.000);2383 
eT 0.09971 +8 0.09971+ 8 0.099717 8 
e4 0.100936 0.100943h); 0.10091,3),0 


(continued on page 19) 
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TABLE 12 (Continued) 9 
OBSERVATION 1 OBSERVATION 2 OBSERVATION 3 


M = (.9016) + (13 .43513)) (t-T> 4+0.0002839(t-To)2 -0.11x1075(t-2, )?-0.000631x1073 


cos W 
(13 43513h) (t-To),, 5381719091 =53.81595623 =53.81323966 
(0.0002839 )(t-T,) 0.00455537 000455516 0.00455)70 
(0.0000011) (t-T,)3 0.00007070 000007070 0200007069 
(0.00000063) cos w 0.00000023 0.00000023 0.00000023 
MT> 029,016 0.9016 0.9016 
M =53 81256461 5328113301,  $~53.80861),0) 
Mt (Revolutions) 0012759539 0.12882986- 061315596 


Mt 45°.93h3040  469.3787)960  7°.35651560 
Mt, 45°56103" $25 46°2201h3"199 79211231. 50) 


4 
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TABLE 13 50 


COMPUTATION OF ECCENTRIC AND TRUE ANOMALIES 


OBSERVATION 1 OBSERVATION 2 OBSERVATION 3 
2 
E=M+esinM+S sin 2M 


2 


e 001018958 0.01018958 0.01018957 
e2/2 0.00509h79 0.00509),79 0.00509h79 
l+e 1.100936 1.1009h3h) 1.100930 
l-e 0.8990565 0.89905656 0.89905660 
(1+e)/(1-e) L 1.225519 1.22u55h1) 1.225504 
(1+e)/(1-e)] 1.10659577 1.1065957h 1.10659570 
sin M 067185327 027239160), 0273558351 
sin 22M 0.999),6819 0.9988h211 0.9966186), 
e sin M 0.0725322h 0.07307h58 0.07425230 
(e2/2)sin 2 M 0.00509208 0.00508889 0.00507756 
E (radians) 0.87932980 0.88762535 0.90585750 
E 5023818863), 5098571863), 51990181159 
E 50°2215",.79. = 50951125", 871 51°95)! 06,522 
E/2 25011127".396 =. 259251)0",986 =. 25°57103",261 
tm E/2 44 047037122 047544692 0 44866723) 
(+e) (1-0) |®tanB/2=tan v/2 0252051080 0.5261275) 0.5385h952 
v/2 27929150".872 27°45'00".338 28918116".789 
V Sh°S9'h1".7hh = 55°30100".676  56°36133".578 
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COMPUTATION OF GEOCENTRIC DECLINATION, GEOCENTRIC RIGHT 
ASCENSION, AND GECCENTRIC DISTANCE OF THE SATELLITE. 


OBSERVATION 1 OBSERVATION 2 OBSERVATION 3 


sin 0 = sin i sin (w+V ) 


(w+v )} 1669222" 6862 166°53'0)".191 16795912" .367 
sin (w+y ) 023550558 0.22691483 0.20799531 
sin 0.51767853 0.54:767853 0.51,767853 
sin 0212898135 0.12),27638 0.11391457 
6} 792h*38" 636 = 7908920".290 6932127" 682 
X= whe! arctan lees i tan (w+y | 
tan (w+y ) =O 6214232136 =0. 2329925), =0 02126591 
cos i 083668886 0.83668886 0.83668886 
cos i tan (w+y) =0 62027):758 =0 619194226 ~0.177918),6 
, 205°52!5".903 205952153".28h 2059521),9" 721 
arctan [cos i tan(w+v ) 168932!19",601 168°58:08".322 169°5)')1".720 
ex Wy°25'1y" 50h =14951'01".606 =15°Y4 7931" dy 
a=q/ (l-e) 
Q (Mm) 6672127351 607212739 607212739 
q (earth's radius) 1205375739 1.05375739 1205375738 
a 


1.17207021 1.17207018 1617207012 


R= a (1 -e cos E) 


e 0.100936 0.10093) 0.1009);3)0 
E 50°22'5h".791  50°51'25".871 51°5)806", 522 
cos E 0,.6376675h 063125552 0.61701099 
e cos E 0.06),36837 0.06372110 0.06228319 
l1-ecosE 0.93563163 0.93627890 0293771681 
R 1.09662596 1.09738h58 1.09906985 
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earth's radius. The value "a" can be computed from the formula: 
asq/(l-e) 

where q is the distamce of the satellite at perigee. The computation 
of these values are shown in Table lh. Cs the Greenwich hour angle 
of the vernal equinox, is tabulated in Table 15. The computations of 
the sidereal times of observations from which this value is derived is 
not shown or listed. The sidereal time was established in the same 
manner as it was for the observations from the known station. 

Now the computed apparent right ascension and declination mst 


be found using the formlae 12, p.37| : 
! 
Xqg sin G + y, cos G 
/ 


9 
ton Oo = 0s (Xe -CX) (R sin O -2,) 
¢ ~ ah X cos G+y, sin G +k ere 


a 


and 


Then the differenres between the computed and the observed values 


are determined by: 


[el = Cee OC 


A6 =6- & 


The values and computations of 5 ¢ and (CX, and the increments 
A a and A 6° are tabulated in Tables 15 and 16. 
The differential equations to be solved [12,5436 | are 3 
aj dxqg+ a2 dya = a3 aX | 
by dx, + be dya+b3 dz, = bj, ak + by a0 
The coefficients of these equations are; 
a, = sin G+ta (1-2) cos G 
cos G = tan (A= CX) sin G 


sec2 (1 -CL)(R cos é - X, cos G + Y,q sin G) 


ac 
mS 





ee —i(< —e- © 


> 
- 
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by = tan 6 cos 6 
Se sin G 
b3 - cos (XK - X) 
b, = sin (X-CX) (R sin O - za) 


/ 
bs = sec® 6 (R eee) “xX, cos G+ ya, sin G) 


be 


The computation of these coefficients are contained in Table 16. The 
solution of the first approximation differential equations is shown in 
Table 17.  ACIC Technical Report No. 86 suggests these equations be 
solved by Doolittle's method but being more familiar with the Gaussian 
solution of normal equations, I decided to use the latter method. The 
solution of the normal and reduced equations yields the increments 

Axa, Aya, ana Aza. These increments are then added to the origin- 
al approximate coordinates to find the corrected position of the observer. 


The corrected space coordinates ares: 


t 06351162753 + 0.912606236 a 0.21092082) 

Since the International Ellipsoid was used as the reference ellip=- 
soid,the xa, ya, Zq coordinates can be miltiplied by 6,378,388 meters, 
the semi-major axis of the ellipsoid, to find the metric distance of P, 
the radius of the observer's position from the earth's center. This has 
not been done because it is of little interest in this problem. How- 
ever, the distance of the observer from the earth's center is found in 
units of the earth's radius by the equation: 

P 2- [x3 + y¥8 428 = 1.000131808 
Then the resultant geographic coordinates of the observer may be easily 
derived as follows: 


sin q’ a ao = + 0.21082979 
tan CD + 0621567762 


SOLUTION FOR conrunn 


Xa 

Ya 

Zs 

Xo sin G 
Ya cos G 
=X. cos G 
Yq sin _G 
R cos 


X,sinG+y, cos G 
-x,cosGtyasinGt+ 
~ y cos / 


tan (X=) 
(At = &) 


Cc 
Cc 


TABLE 15 


5 
i 
APPARENT RIGHT ASCENSION, s 


APPARENT DECLINATIQN, 


FIRST APPROXIMATION 


OBSERVATION 1 


14,°2511)" 50h 
792h! 38.636 
1.09662596 
355° 6t30". 
=30° 51100". 
85~31°54".975 


=@71906' 0" 671 
~0 91,614,889 
Osseo 13170. 
0012898135 
0.9916),703 
32922915" 1:96 
031971509 
0294751372 

—0 3372529 
1.1138558), 

=0 659730303 
1.35677091 


tan (CX, -Q{) = 


=X, cos G+ ya sin G+R cos } 





—_——___-_> 
ee 


-0.33303764 
0.29560113 

~ 011395127 
=0 86393336 
1.087),6588 


=0 0373651 
0210958125 


=0.0341632h) 
342908117" 06 
355°33*31" 910 


/ 
tan 6,5 


OBSERVATION 2 


1°51? 01' .606 
7°08120" 290 
1.0973858 
359° 391% OO". 
=32° 521 OO", 
85°33 *5kt 6405 


=70°,2 8 Se" ° 799 
=0 94388553 
0.33027279 
0.12)27638 
0.9922h765 
34°47? 58" 39k 
~0 »26219670 
0.96501)95 
=0..27170236 
1.07382217 
~0.6),610,06 
10417450),6 


OBSERVATION 3 


1597931" hyd 
693227" 682 
1.09906985 
8~ 03* OO". 
=36° 39! OO", 
85°38 17" .160 


69950 #45" .719 

~0.938770L, 
0 .34h5h),07 
011391157 
09939055 

= 7h? 31" oy 
=0 137137) 
0.99088)56 
~0.13595301 
1.0181,8322 

=-0.7)4),020),0 
1255356636 


x, sin G+y, cos G © 


+ 0235199282 
+ 0.913105082 
+ 0.210569117 


=0.3322095 
030157376 
011625366 
~0 86186667 
1.08887727 


=() ge 030667 19 
0.1107569) ~ 


= » 27688730 
34°31! 23" 0757 
359922825" 363 


— 


<—______—_ 
ee 


~0.330)),037 

0.311609, 
=—0 12127705 
—0 85719579 


1.09191551 


=~0 01583543 
0.113))267 


=~0 13958972 
=7 9567" 593 
7950°h3" 8448 


cos Coe ~Q) (R sin O - Za) 


=X, cos G+y, sin G+R cos 





TABLE 15 (Continued) 


FIRST APPROXIMATION 


OBSERVATION 1 
/ 
cos (CX .- &) 0.94630093 
R sin 0.141430 
(R sin () -z,) -0.06912),82 
“X, cos G+y, sin Gt 
R cos 0.10958125 
cos( @! 


7X) 
(R sin O-za)  -0.065),1288 


tan Og ~0.59693197 
Cc 


30°50! 0h" .037 


OBSERVATION 2 


0.96373889 
0.13637898 
0.1107569), 


-0.07149992 


-—0.64555702 


=-32°950!)0" .376 


OBSERVATION 3 


0.9903976 
0.12520007 
~0.08536905 
0.1134)267 


=0 085929 


-0,74530)13 


~36°)1150" 335 


De 





TABLE 16 


COMPUTATION OF COEFFICIENTS OF DIFFERENTIAL EQUATIONS 


A&G (radians) 
sin G 7 
tan( CZ -COcos G 
al 


cos G ' 
tan(C¥-CQ)sin G 
a2 


sec?( (X'-CX) 

R cos 

X, COS G 

Ya sin G 

Rcos A-xa cos G+ 
Ya sin G 


a3 
Pl 
be 
B3 


sin (CX'-C) 
R sin 


Za 
R sin -Z 


FIRST APPROXIMATION 


OBSERVATION 1 OBSERVATION 2 OBSERVATION 3 


ay = sin G+tan (= (CX) cos G 


0°13158" .090 0°1613h" 637 0°12916".152 

0° ,.23280278 0° .27628805 0° . 20h1:8666 

0.00L,06318 0.008221) 0.00356897 
~0°O1? 0)" .037 ~0°01!19" 62h 0002150" 335 
~0° 01778806 ~0° .02211778 0° 0731527 
-~0,000310))6 -0.00038603 0.00082581 
~0.9461889 094388553 ~0.9387701 
—0.10923529 -0.08973590 -0,0168)180 
~1.05538)18 ~1.033621))3 -0.9856119); 


a9 = cos G - tan (- CL) sinc 


0.32373178 0.33027279 0.3): 5407 
0.31925)56 0.2561:5593 0.12762863 
0.00hh7722 0.07381686 0.216915) 


ay = sece ('-CO(R cos()=Xa cos G+y, sin G 


1.1138558) 1.07382217 1.018)8322 
1.087),6588 1.08887727 1,09191551 
0.11395127 0.11625366 0.12127705 
-0 .86393336 -0 86186667 -0.85717579 
0210958125 01107569); 0.113)267 
012205772 011893326 021155396 


bz = tan 6! cos G 


-0.19336597 ~0 21339059 ~0 625634782 
{ 
bo == mA sin G 
~0 56513760 ~0 6098827 ~0.698)6h).1y 
b3 = -— cos (CL GZ) 
~0 9751372 ~0.96501L45 =0.99088)56 
bj, = sin (O'-a )(R ao ~Za) 
~0.31971509 ~0 26219670 ~0.137137h 
Oo UWydhh30 0.13637898 012520007 
0.21056912 0.21056912 0.21056912 
~0,06912))82 -0.071901y =0.08536905 
0.02210025 0.019452h1 0.01150038 


oT 
TABLE 16 (Continued) 


FIRST APPROXIMATION 


OBSERVATION 1 OBSERVATION 2 OBSERVATION 3 


f 
by = sec* 6 (R cos b ~ xX, cos G+y, sin G) 


} 
sec2 1635677091 1641745046 1255356636 
R cos - X, COs G+ 

Yq sin G 0.10958125 0.1107569i 021131267 


bg 021867665 0.156992),8 0.1762),072 





— eee 


— 


dxa 
dya 
az 


TABIE 17 
SOLUTION OF DIFFERENTIAL EQUATIONS 


FIRST APPROXIMATION 


58 


Yo+t 02912606236 


OBSERVATION 1 OBSERVATION 2 OBSERVATION 3 
f 
ax AQ 0.0004959L 0.00057351 0.0001236 
bj, AC 0.00008980 0.00009380 0.000010) 
bs A 6! -0.0000616 -0.00006060 0.0001455) 
bj, Ac +s A 6' 0.000036), -0.00003320 0.00018658 
OBSERVATION EQUATIONS 
ajdxq + apdy, = a; dCr’ | 
bydx, + body, + b3 dZ,q = by, aC! + bod 
dXo dyag dZ. L = 
~1.05538)18 0.007722 ~0.000959) =-1.05140290 
~0.19336597 ~-0.56513760 -0.914751372  -0.000036) =1.70606093 
-1.033621))3 0.07381686 ~0.00057351 -0.96037808 
~0.21339059 -0.60984827  ~0.96501h5 -0.00003320 -1.78828651 
-0.9856119)) 0.216915); =~0.0001236 -0.76910886 
~0.2563:782  -0.6986h14  -0.99088456 -0.00018658  -1.9):588310 
NORMAL EQUATIONS 
dxa dy¥a dz L y 
3302280071 0012365699  0.643153009 0.001585979 4.07066)758 
———s 1.231668857 1.816085355  0.000041227 3.171h1137 
22810887350 0.000258267 5.270383981 
~~ §,000000783 0.001886255 
REDUCED EQUATIONS © 
dx, qdYa dz g £ ys 
32302280071 0.12365699  0.643153009 0.001585979 .07066)758 
, 3  1.227039250 1.79200)08) ~-0.000018156 3.019025178 
-0°19):760285 ~1.460429309  0.068531)01 -0.000024103 0.068507298 
-0.000):80268 0.0000311797 0.00 707 02000000013 0.000000012 
dXq dy dZa 
~0.000530089 -0.000f98816 +0.000351707 
x 0.3519928)2 V 0.913105082 Z 0.210569117 
Axa -0.000530089 Aya ~0.0001988)6 Aza 0000351707 


Zo+ 0021092082) 





59 
/ 
tan o= tan Daze? = 0.21567762/0.99327733 = 021713736 
tan A = Yo/x. = 20596594K6 


P 
A 


The differential method was repeated again for a second approxima= 


which result ing 
12° 15% 03".168 N 
68° 56% 1)",336 W 


{| 


tion. The corrected space coordinates from the first differential 
process were used as the approximate observer's position. The data 
used in and the computations of the second differential approximation 
are listed in Tables 18, 19, and 20. 


The corrected and final results of the second approximation are: 


+ 035188350) aa 0.913709058 + 0, 290819193 
Pp = x6 a Yo° + Zo* = 1.000517795 
sin p! = 20/) = 0.2107),007 
tan ('= 0.21558157 
tan D = tan P/i-e2 = 0.21558157/0.99327733 = 0621704066 
tan = Yo/x, = 2059673370 
The corrected and final geographic coordinates are: 
D = 12° 1h! h".118 W 
A = 68° 56! 18".0h5 Ww 
Further approximations would probably yield changes to these 
coordinates but as stated previously the assumption that two approxi-~ 
mations are sufficient has been made. 
As a note of interest, and to satisfy my own curiosity, further 
calculations were made using the linear method but are not included in 


this paper. First, I used my own computed value of the Gaussian gravi- 


— 


tational constant, k, in the formula to determine "a": 








—- ee 


TABLE 18 60 


(APPARENT RIGHT ASCENSION, Xe 
SECOND APPROXIMATION 


OBSERVATION lL OBSERVATION 2 OBSERVATION 3 


tan (Qi-) _ Xa sin G + ye cos G 
~X,g cosG+y, sinG+R cos 
Xa ma + 0.351462753 —_— 
Ja Se + 0912606236 —— 
Za ina + O¢ 21092082) a 
X,sin G =0 633253609 =0 3317060 —003299427h 
yacos G 002954396) 0.3010901 0031443307 
=—Xo cos G ~0.11377966 ~0.11607858 ~0,12109h))3 
Yq sin G ~0 86316138 086139582 ~0.856727h9 
R cos 4 1.0876588 1.08887727 1.09191551 
XqsinG +y,cosG 0 .03709645 ~0.03033159 -0.01550967 
-XgcoOsG+y_ sin G+ 
R cos 0.110228) 0.111),0287 0.1109362 
tan (OX -@& -0.3365525), -0. 27226938 ~0.13593811 
(COO 341923°57".123 3h 6409".492 9 =7°W 228" 2) 
Ge 355°H9911" 627 = 359937111".098 8° 03#03".017 


tan on 2 cos(CX.-(X )(R sin 6 Za) 


- Xg Cos G+yYag sin G+R cos 


cos (X-.-CX) 0294776396 0.96)87589 0.99088653 
R sin 05141430 013637898 012520007 
R sin§ =z 5 =0 .069),7652 ~0 074518) -0 208572075 
cos(Q'=O9(R sinOnzq) -0.065873 =0.07192362  -0,0893951 
tan b% —0.59739111 =0 6615561730 =) 6 7))}7237 


4 
OF =30°51113" 9391 


~32950'h9" .152 


=36939! 59".995 





TABLE 19 61 
COMPUTATION OF COEFFICIENTS OF DIFFERENTIAL EQUATIONS 
SECOND APPROXIMATION 
OBSERVATION 1 OBSERVATION 2 OBSERVATION 3 


a, =sin G + tan Coane ) cos G 


CeO NO! ~0°0211",627  0°01')8".912 -0°00!03".017 
CX! =0°.0448)639  0°.03025333 + =09,00083806 
| (radians) =0 00078272 0.00052802 ~0 00001463 
/ f_ 
6 -6,=AQ 0°00913",391  ~0901'10".8h8 0°00159".995 
0°.00371972 ~0°.01968000  0°.01666528 
Aé (radians) 0.00006),92 ~0 0003381 000029086 
al -1,05538)18 ~1,03362143 ~0.98561194, 
a2 = cosG-tan (C'-® sing 
a, 0047722 0.07381686 002169154) 
a3 = secé (CV'-X)(R cos 6 - X, cos Gt+tya sinG 
sec’ (X'-CY) 1.1138558), 1,07382217 1.018),8322 
R cos) 1.0874,6588 1.08887727 1.09191551 
X_ cos G 0211377966 0.11607858 0012109441 
y., sinG =0,863),6138 =0 86139582 0 .856727h9 
(R cos b=Xa cos G+ 
Ya sin G) 0.110228), 0.111),0287 0011409361 
ay 0.12277h58 0211962687 0.116202),3 
by Sea cos G 
by ~0.19336597 =0.21339059 ~0.2563),782 
bg =-tan 6 Sin G 
bo ~0 56513760 0260984827  -0.698)6)1) 
b3 = = cos (CU =O 
b3 =0 94751372 =0 9650145 ~0.99088)56 
bj, = sin (Cx'-CX) (R sind ~%.) 
sin (X'-%) =031971509 =0026219670  =021347137h 
R sin Oo Lyihhh30 0.13637898 0212520007 
Zo 0.21092082 0.21092082 0.21092082 
R sin () = za =0 0697652 ~0.0745418 =0.08572075 


bh 0202221269 00019562 0.0115h776 


TABLE 19 (Continued) 62 
OBSERVATION 1 OBSERVATION 2 OBSERVATION 3 
bg = sec® 5'(R cos ()-xa cos G+yqg sin G) 
seo? 135677091 LoL 7HSOK6 —=»-1459356636 


(R cos 6 -X_ cos G+ 


Ya sin G) 0.110228) 0.111,0287 0.11),09361 
be 0.1495h986 0.15790805 017725199 





TABLE 20 63 
SOLUTION OF DIFFERENTIAL EQUATIONS 


SECOND APPROXIMATION 


OBSERVATION 1 OBSERVATION 2 OBSERVATION 3 
ea =0.00009610 0.00006317 -0,00000170 
vpA -0.00001739 0.00001032 =-Q.00000170 
Beko! 0.00000971  =0,00005),29 0.00005155 
bP AMX! + be A 6 =( 600000768 =-0.0000);397 0.00005138 
OBSERVATION EQUATIONS 
al Xa + a2 dyg = a, dC! ' | 
by dxg + b2 dy, + b3 dzg =b, dX + be ad 
dxe dys dz, 4 s 
~1.05538),18 0.00L7722 0.00009610 -1.05081086 
-0.19336597 =-0.56513760  -0.9h751372 0.00000768 -1.70600961 
~1.033621)3 0.07381686 ~0.00006317 =0.9598677) 
-0.21339059 +~0.6098827  =0.96501h5 0.0000)397 =-1.7882093) 
~0.9856119) 0.216915); 0.00000170 -0.76869),80 
—0.25634782  =0.698)6h1i -0.99088456 -0.00005138 ~1.94571790 
NORMAL EQUATIONS 
dxa dya azo £ th 
3302280071 0.123645699  0.643153009 ~0.000035501 .06903279 
Seen. Me. 1.231668857 1.816085355 0.000000868  3.17100778 
: ~ 2810887350 0.000001203 5.270126916 
oo  -—-Q,.000000018 -0.000033)12 
REDUCED EQUATIONS 
dX qdya 0% 5 Me D. 
xg 3.302280071  0.123645699  0.63153009 -0.000035501  .06903279 
dyg -0.037442523  1.227039250 1.79200h084 0.000002197  3.0190)5531 
AZe -0.19:760285 -1.460429309 0.068531)01  0.00000L909  0.068536309 
0.000010750 -0.000001790 ~-0.000071603L | 0.000000017 0.000000016 
dXa dya dz. 
+0.000020851 + 0.000102822 -0.000071631 
Xg 0.351)62753 Ya _0.912606236 Za 0021092082) 
A\Xat 0000020851 Ya + 0.000102822 /[\Za ~02.000071631 
Xo + 02351483504 Yo + 00912709058 Zo + 0.2108),9193 





él 
Te’ 


n 


where n is the velocity of the mean anomaly, M. The resultant "a" from 
this formula has the value 1.17190297 in contrast with the value 
1.17207018 which was computed from the formula a =aq/l-e. The final 
result of geographic coordinates using this value and the linear method 
is listed in the summary. 

Second, the linear method was used again in which the eccentricity 
derived from the mean orbital elements was used instead of the smooth 
value for e. The mean eccentricity for the time of observation is 
0.10067497. The resultant coordinates from this experiment are also 


listed in the summary. 





Summary: 

Method Latitude , PD Longitude , A 
Approximate 12° 13! 23".21 N 68° 55* 07".86 
Differential (1st) 12° 15° 03".168 68° 56% 14".336 
Differential (2nd) 12° 1h2 bh" ,118 68° 56* 18.045 


Linear (using computed k) 12° 15% 06",351 68° 560 11".075 
Linear (with mean e) 12° 15° 15" ,276 68° 57: 25.278 
Actual 12° 05? son.) 68° 508 1h".0 


— aawe €8 f 


_—«<- —_ a (ee 





6. THE INTER-CONTINENTAL TIE 
The last step of this thesis is to effect the "tie" between the 
two stations, the geographic coordinates of which are now known. This 
problem can be considered to be one of geometric geodesy or what is 
more commonly known as the "determination of long lines" and is considered 
to be a part of higher geodesy. 
H. Fo Rainsford states there are two main problems in the field of 
long lines 6, pel2|s 
a) Given the geographical coordinates of two points, 
to find the distance and azimuths between them; 
commonly called the reverse (or inverse) problem; 
b) Given the geographical coordinates of one point 
and the distance and azimuth from it to another, 
to find the coordinates of the second point; 
commonly called the direct problem. 
Since only the geographic coordinates of the two points are known 


in this case, the reverse problem applies. 


phe following data are givens: 


Point 1 Point 2 
aC a0, es Olifantsfontein, South Africa 
: “®, 12° Ly! bht.118 N dh, 25° 57° 3h'.700 S 
—— X, 68° 568 18".0h5 W A, 28° 1s! 51".100 E 


Problem: Find the forward azimth, Xj, at Point 1; find the 
back azimuth, 915 at Point 23 and determine the distance, s, between 
the two points. 

Since the reference surface chosen for this paper is the Inter- 
national Ellipsoid, the distance, s, between the points will be one of 
the basic geodetic curves called the geodesic. It is defined as l13,p.2]¢ 

a line on any mathematically defined geodetic reference 

surface that defines the shortest distance between amy 

two points on that reference surface. It is. geometrically 

equivalent to a great circle on the sphere. It is a wmique 

line on the reference ellipsoid that contains the normal 

to the ellipsoid at each point. 
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66 
Refer to Figure 7. It canbe seen that the line, s, connecting the 
two points, Pl and Py, contains the normal to the ellipsoid along its 
entirety. The forward and back azimths, 12 and X97 respectively, 
are shorm at the points. The curvature of the geodesic in the diagram 


is exaggerated for clarity of figure only. 





Figure To 


It should be noted that the distance difference between a geodesic 
and a plane curve is very small but the difference in the azimuths of 
the two would be appreciable. As the distance between points increases, 
the difference between the two lines (geodesic and plane curve) diverges 
rapidly. 


There are several methods available in geometric geodesy to deter= 


57 
mine the geodesic and the azimuths but because this particular problem 
is concerned with the inverse problem and the distance between the two 
points is approximately 6,000 miles, then a choice of only two methods 
is available. These are the Rainsford Method and Sodano's Fourth 
Method. Another method that solves for distances greater than 5,000 
miles is Rudoe's Method Pie Since this procedure treats the normal 
section only, it is not applicable for this problem. The method 
chosen was Sodano's Fourth which is a non-iterative solution. The 
accuracy is of the order of +0".1 in azimuth and + 0.5 meters at 


P 





Sphere 


Figure 8. 
Proceeding with the solution, an auxiliary spherical triangle is 
employed. In our problem, the westward point, P], is Curacao and the 
eastward point, Po is Olifantsfontein, South Africa. The difference in 
longitude, L, is always positive. /33, poll | ° 
The first step is to find the parametric latitudes, Bi ana B 2. 


using the formula 


tan [3 = (1 =f) cee 
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with the results 
Bri = 12° 078 33.522 
B 2 == 259 53" ODM L8e 
Using the parametric latitudes, the approximate spherical distance, 0 ’ 
is found from [13,p.la | 2 
cos G = sin fi sin [Ep + cos B, sea), cos L ==0.2017305506 


and 
2 
*= 0.77515)8582 


sin Gi = |cos fi ie B, sin’ $ + sin’ A A 


The resulting value for G is 101° 38! 17".336. 


Now X must be found using the formula ee L+( , = L) where [13 ,p2h2] : 


wr -»=|6 ap B 0 at: orn sin 6 + (F+E) G |e" aid 0°18!20" 866 


2 
where DP = 20626)''.80625. The value for Me 
97° 11! 09" .1h5 


0° 18% 20" .866 
97° 29! 30".011 


WoW 


L 
(Ee) 
r 
Then the forward azimuth {19 is computed from the formula 


(23, pe2| 3 
cot 49> [tan 8, con), - cos i sin i] fein = 12,°16!06" 607 
The back azimuth is [ 43, peli2 | 2 
cot X91 [sin 2, cca = cos [3y, tan &, | fos A. = 277°)9356" 50h 
Sueincatisten ie to find the true spherical distance, Gq usiiewe— 
97°29130".011 instead of L = 97°11'09".145 in the same formula as 
before: 

cos 6. = sinfJ, sin [3, + cos B walls cos A = -0.2063866912 

GO, = 101° Shr 38".472 


The final step is to find the geodesic distance, s, from the formla 


[13, p.l3 | 3 
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s =b(Ao Q,+ Bo sin = cos 20 - Co sin 2 oe cos h GO ) 
s = 11,309,342.6228 meters = 6,102.21 nautical miles. 
Refer to the following tables for the detailed computations of the 
constants, coefficients, and the unknowns. 


Summarizing the results: 


Point #1 Point #2 
ers NW... Olifantsfontein, South Africa 
Ob .118 N an 25° 57% 3".700 S 
i, 68° 56118" O45 W X, 28° 1h 51".100 E 
Qio = 114° 16* 06".607 
CL ai = 277° ko" 56".50h 


Ss 11,309,3h2.6228 meters. 


This concludes the computations. 


TABLE 21 70 
LIST OF CONSTANTS 


FOR USE WITH SODANO'S FOURTH METHOD 


INTERNATIONAL ELLIPSOID 


/Q= 297 00 
a = 6,378,388.0000 meters - 
b = 6,356,911.9)62 meters | 
e = 0.0819918898 e¢ = 0.0067226700 
dd = 0.0822688896 e! 2 = 0.0067681702 
i 
7 © 
e’ + e 


| 2 
Kj = ed = 237.2388918 


i yee 16 = 7O519.51 15 
7 2(16e2n¢ +e 2 : 


e€ 


Ky = 1e2 Wete *  _ 1, o86s206h9 


el 
aes “ 
lée2 Ne 
Ke = = 007994593686 
ome 
¢ 16e°N2 + & 


Kg = 1822 N2___ _ 5, 986820609 
e 


TABLE 22 
COMPUTATION OF VARIOUS TERMS 


ea Jas, cos i sin L 


sin 4 
cosh cos Bu 08796086600 

in 029921455600 
ion cos 0.8726998266 
(cos iS cos Best sin “ty sin § 0.8910182190 

B = A& = 0.793913h666 

. = (cos ~ B cos 0 i 
K3 b,.9865206))9 
B cos ~0.1601566007 
(cos = B cos @ ) =0..0415739h99 
(cos 9 -B cos 9 )K3 = 0.008337 2662 

D= -K, sin J, sin), 
0.1:010812630 
stn 8, sin By, -0.0917015172 
—K), sin , sin}, = 0.0367797603 
BE = -Kg sin/}; sin Bo. 

Ke 5 ein eee 
Ss By sin By =0 0917015172 
-Kg sin Bi sin}, : 0 .0733116370 

. F = CK6 
C 0 .0083372662 
K6 398652069 
CK6 = ~0 .0332366839 
— fossa 
G = 9°/sin§ 

6 9 101938117" .336 

, (radians) 1.7739203h4 
0 3 616793387 


6°/ sin 0 = 3,.2128)6105 


fo 





TABLE 23 


COMPUTATION 


Y ( 


PARAMETRIC LATITUDE & APPROXIMATE SPHERICAL DISTANCE 


Parametric Latitude ( B) 


tan 9 = (1-£) tan 


0.003367003), 
(1-f) 09966329966 
tan D 021558157 
tan B 02148557061 


f=12° 07133" 521 


APPROXIMATE SPHERICAL DISTANCE 


— 
—_ 


—0 .1:8686087 
~0.852216078 
-25° 53101" ,189 


cos 6 = = sin 8 sin, + cos B, cos J, cos L 


=0 43654528h)9 
0.8996822851 


tr =A 979 11109" hs 

sin L 0099214556 

cos L ~0 12508862 

sin 0.2100618662 

cos 0.9776880956 

me ow —Q o0917 015172 

a | cos 77 0.8796086600 

cos /Y cos B cos ‘L =0 110029033), 

cos @ ~0.2017305506 Be: 2 

in G [cos {3, cos [3 sin 3 + sin® | Le At) * 

h8 eaten 
0.75002937 
0.5625)))0559 


in 2 
. Lp =38°00! 3;"".710 
=19°00! 17" 0355 


R 0 63256771 
ae “RY 0.1060)6),310 
cos [j, cos /§2 sin@ L/2 0019)48186232 
eer ug oe 4+ sin? 
0.60086505)2 
sin 0/2 067751548582 
A/2 50°h9'08" .668 
O0 <= 10193817" .336 
sin Q) 0097941057) 
cos a =() » 2017305506 


TABLE 2h © 
COMPUTATION OF ( = 1), Kyo and X 9). 


(et) = (Kk) + B ) 6 + (¢+D) sin oo ae Ap 


Ke 

Ri 237.2388918 
K2 70519 51145 
(Ki+B) 2380328052666 
(Ki +B)O 4.22..25123586 
(° +} 0.028429) 
C+.D) sin8 0.0278h774.65 
(F + EB) 0.0)007),9531 
(F+E) G 0.128751,6570 


[ j 
Qen eee ndilieade 0": 2 Ea 


A=L4+( Je L) = 97° 121091, 1,5+ 0018? 20".866 = 97° 29°30".,011 
cot Xio = ltan B,cos Be cos A sin | / sin A 


cot X40_ (-0.1)852216078 )( .9776880956 )—(-.1303820)28) (0.2100618662 ] 
7 0.991,638331 


cot Kyo9= =—0.)70070945/0.9914638331 = -0.))508556737 
(ip = 12)° 16 06".607 


cot 4 = [sin [3, cos A - cos £, tan B,| / sin Ne 
cot Xo) — | (-.43651)52819) (-0.1303820)28 )~( .8996822851) ( .2148557061) 
—= | 00lti— ——_—— 
cotKoz= -0.13638)2066/0.991)638331 = - 0.137558)282 
(X91 277° 9 56".50) 





TABLE 25 7h 
COMPUTATION OF TRUE SPHERICAL DISTANCE, 0. 
(use Avice 1) 


COMPUTATION OF GEODESIC 


cos @, = sin / sin 4+ cos #4; cos #2c0s A 
= (=.0917015172) + ( .8796086600 )(-.1303820);28 ) 
= (=.0917015172) =(.1168517)0) 
= =,2063866912 
GO, = 101° Skt 38 
Oo = 1.778677025 radians 
cos (a = [cos | COs By, sin | /sin Q, 
cos 5) cos Br 0.8796086600 
sin‘) 09914638331 
sin 6, 0.978)705080 
cos 0.8912891769 = maximum latitude 
Bo of geodesic 


bo =1+ (8 %/y) sin? QZ, - (384) sink, 
é _. po mn 
0 62056036030 
sé 00000021473 
sin 0.0)422728h16 
Ao 1.000347 7993 
~ (© ®/u) sin? hi, -(é 4/16) sint 
eh 00016920426 
2 Po 02026036030 
Mag 00000028630 
ink RB, 0.0227 28h16 
“A. 0.0003);77691 
/ 
e hk 
@ {328 00000003579 
sink oe 0.022728)16 
Co 6.0000000151 
cos 2 O —_ 2 sin oF sin 2, _—- cos 6 - 
sin@/} 
2 gen 0.4.20123732h 
fp 2 =0 9136545289 
aie G ee 02056036030 
cos Bo 02063866912 
cos 2 =0 .6856357822 
2 133°17 #08" 261 
ho 26603), 16" ,522 


cos hO 


cof) ° 059807 1587 
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TABLE 25 (Continued) 


s = b (Ap 0, + Bo sin GO, cos 20 ~ Co sin 2 O, cos h O ) 


Ao 1.0003h77993 
6 o 12778677025 
A, Oo 107792956471 
Bo 0.0003477691 
sin 0, 0978705080 
cos 20 ~0.6856357822 
B sin), cos 20 0000233309) 
BA 0.0000000151 
sin 2 Oo ~0 .1:03886581); 
cos }; 00598071587 
Cc, sin 2 6, cos hO 0.000000000 
b 6,356, 911.9462 
s 11, 309, 32.6228 meters 


6,102 .hh21 nautical miles 
Py02T woes statute miles 


To CONCLUSION 

Upon reviewing the results of the determination of the orbit of 
Explorer I, it can be stated that, in general, reliable orbit para= 
meters ‘cannot be reduced from only one set of three observations of 
the satellite. 

Based on the results of the calculations of i ans, the position 
of the orbital plane in space can be accurately determined from a few 
observations. To determine the size and shape of the elliptical orbit 
and to accurately orient the orbit in its plane requires a series of 
precisely observed and reduced coordinates. As Moulton states |b,pol] : 
the difficulties arise in finding a, e, and w. As mentioned previous= 
ly, these three parameters determine the size, shape, and orientation 
of the orbit in its plane. To determine these three parameters satis~ 
factorily, the whole time interval of the three observations should 
be divided into two nearly equal parts by the second observation. If 
this is so then 01 and (3 will be nearly equal. With the time inter- 
vals equal or nearly so, the ratio of the area of the sector to the 
area of the triangles (see Figure 5) will be more reliable and the 
series which is used to determine the auxiliary parameter p will 
converge more rapidly. Additional terms can be used in the series 
to determine p. 

In an analysis of the solution for ul, ug, and u3 (derived from 
Pp), upon which 0. 9 oe 959 and O 3 are based, the reader should note the 


value of each u determined by its sine and cosines 


Uy ud U3 
gine value 248° 13102".933 253. 13?50".317 255° 3,925,033 
cosine value 21,80 13°02" 756 253 13250" ,305 2550 ),6:17" 81h 
difference 0" .177 0" ,OL2 11°52" 761 
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us 
There is a close agreement between the value of the angle determined 
by the sine and that determined by the cosine for uz and uy but there 
is an 11'52",781 disagreement between the value for v3 derived from 
its sine and cosine. The value used for u3 in the ensuing calcula= 
tions was an average of these two which is 255° 0?21".h2h. This could 
be the source of the error existing between the computed a, e, and w 
and those calculated from the 82 observations. I believe that this 
can be attributed to the uneven time interval of the observations: 
h 

PSUS} mar.ae 

T3 22 32 03 .100 ) 058.680 
and the unrealistic observed declinations: 


Obs o #1 Obs #2 Obs. #3 
7° 22! S 51° 06: § 9° 25' s 


It seems logical to assume that the declination should either increase 
or decrease progressively during the short intervals of the observations, 
particularly during the same transit of the satellite. There is no 
such progression in this caseo 

The assumption that the orbit is correct is another source of 
error. The accuracy of the orbital elements predicted for the time 
of the observations from the unknown station depend on the accuracy of 
the initial values of the parameters determined from the 82 observations. 
The standard error of one observation of + 5.3 undoubtedly lessens 
the reliability of these values, especially the mean anomaly, M, which 
is probably the crux of the entire problem. 

The assumption that the observed coordinates of the satellite 
and the computed orbital elements are referred to the same mean coordi- 


nate system also influences the result. As stated in | 9op.2 | on the 
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orbital elements: 

the reference plane is defined as the true equator of 

the date. The origin of right ascension is a line 

shifted from the mean equinox of the date by an amount 

equal to the precession in right ascension between 

1950.0 and the date. 

Given below is a formula with which values can be obtained to correct 
the right ascension given in the orbital elements, in a right ascension 
referring to the mean equinox of the date: 

° (7)= (° (or) +3°.508 x 10-5 (Smithsonian Day - 33286) 
where DOI indicates the values determined by the Differential Orbit 
Improvement Program. This information was not available. 

It is apparent that the methods and procedures demonstrated 
herein to effect the tie between the taro points on the surface of the 
earth can, with refinements, be employed in practical applications. 
The success of this application depends substantially on the charac= 
teristics of the satellite itself. In fact, the only satellite that 
would be practical for precise results would be a geodetic satellite = 
one that would have the following characteristics: 

a Have a nearly circular orbit or one that has a 

very small ellipticity -— this would permit reliable 
predictions to be made in its position and eventually 


an ephemeris could be made available; 


be Be spherical in shape - so that its motion in space 
would not affect its presentation surface; 


co Have a perigee of approximately 1,000 miles to 
minimize the effect of atmospheric drag; 


do Have a large inclination angle so that it would 
be visible from most latitudes; 


eo Have a flashing beacon that would be visible at 
its long range. 


Wr 


pak =a ——— = 
On 
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